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[i]  Under  steady  wind  forcing,  wave  development  follows  the  duration-  and  fetch-limited 
growth  laws.  These  growth  functions  are  used  extensively  to  obtain  the  sea  state 
information  when  only  limited  observations  of  the  environmental  variables  are  available. 

Validation  and  verification  of  wave  models  also  employ  numerical  experiments  of 
duration-  and  fetch-limited  wave  growth  as  benchmark  tests.  The  reference  wind  speed 
reported  in  most  of  the  wave-growth  data  is  the  equivalent  neutral  wind  speed  at  1 0  m 
elevation,  Ul0.  It  is  generally  believed  that  a  more  suitable  scaling  wind  speed  is  either  the 
wind  friction  velocity,  u *,  or  the  wind  speed  at  an  elevation  proportional  to  the 
wavelength  of  the  ocean  surface  fluctuation,  U-^ 2-  The  connection  among  the  growth 
functions  using  different  velocity  scales  is  the  drag  coefficient  of  the  ocean  surface.  In  this 
paper,  the  similarity  relation  of  the  drag  coefficient  based  on  wavelength  scaling  is 
applied  to  the  conversion  of  the  wave  growth  functions  from  t/10  to  U\j2  and  w*  scaling. 

The  results  are  in  good  agreement  with  field  measurements  that  include  direct  w* 
measurements.  Comparisons  with  numerical  model  output  are  also  described. 

Citation:  Hwang,  P.  A.  (2006),  Duration-  and  fetch-limited  growth  functions  of  wind-generated  waves  parameterized  with  three 
different  scaling  wind  velocities, ./  Geophys.  Res.,  Ill,  C02005,  doi;10.1029/2005JC003180. 


1.  Introduction 

[2]  Duration-  and  fetch-limited  wave-growth  functions 
are  of  fundamental  interest  in  basic  research  and  engineer¬ 
ing  applications.  They  quantify  the  wave  evolution  under 
various  driving  forces  represented  by  the  source  terms  of  the 
action  density  conservation  equation.  Some  of  the  source 
functions  are  still  poorly  resolved  and  one  of  the  methods  to 
gauge  the  performance  of  numerical  wave  models  is  to 
compare  the  model  results  with  fetch-  or  duration-limited 
wave-growth  functions  [e.g.,  Komen  et  al,  1984,  1994; 
Janssen  et  al.,  1994;  Janssen,  2004],  In  the  study  of  ocean 
surface  wind  stress,  the  growth  function  expressed  as  the 
dimensionless  dependence  of  wave  variance  on  wave 
frequency  is  invoked  to  make  comparison  among  different 
expressions  of  the  drag  coefficient  or  dynamic  roughness 
[e.g.,  Toba  et  al.,  2001;  Jones  and  Toba,  2001;  Hwang, 
2005b],  Wave  growth  data  are  usually  reported  with  the 
neutral  wind  speed  at  10  m  elevation,  {/]0,  serving  as  the 
scaling  wind  velocity.  So  far,  there  have  been  many  differ¬ 
ent  growth  functions  proposed.  The  discrepancies  among 
the  proposed  functions  are  caused  mainly  by  the  stability 
conditions,  the  combination  of  field  and  laboratory  data  in 
some  of  the  analyses  (e.g.,  see  reviews  by  Kahma  and 
Calkoen  [1992,  1994]  and  Young  [1999]),  the  spatial  vari- 
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ability  of  the  wind  field  caused  by  the  proximity  of  land 
[Dobson  et  al.,  1989;  Donelan  et  al.,  1992],  and  the 
difference  in  wave  development  stages  in  individual  data 
sets  [Hwang  and  Wang,  2004].  After  sorting  out  the  stability 
conditions,  excluding  laboratory  data  from  the  analysis,  and 
using  the  average  wind  speed  between  measuring  stations  as 
the  scaling  wind  velocity,  Kahma  and  Calkoen  [1992,  1994] 
found  that  many  of  the  discrepancies  can  be  reconciled. 
Three  sets  of  equations  were  presented  for  stable  stratifica¬ 
tion,  unstable  stratification,  and  composite  data  set. 

[3]  During  the  course  of  growth  function  investigation,  an 
interesting  question  was  raised  regarding  the  best  scaling 
wind  velocity.  The  adaptation  of  C/jo  as  the  reference  wind 
speed  is  mainly  based  on  practical  considerations.  The 
scaling  laws  call  for  “free-stream”  velocity,  Ux  [e.g., 
Schlichting,  1968;  Kitaigorodskii,  1973].  For  the  marine 
atmospheric  boundary  layer  modified  by  the  ocean  surface 
waves,  one  expects  that  the  dynamic  influence  of  surface 
waves  to  decay  exponentially  away  from  the  air-sea  inter¬ 
face  with  the  decay  rate  scaled  by  the  characteristic  wave¬ 
length  of  the  surface  fluctuation  [e.g..  Miles,  1957;  Phillips, 
1977].  It  is  logical  to  consider  that  the  wind  speed  at  an 
elevation  equal  to  one-half  of  the  characteristic  wavelength, 
Ux/2,  to  be  a  reasonable  representation  of  U0 0.  Another  good 
candidate  for  the  reference  wind  speed  is  the  wind  friction 
velocity,  tr*,  which  is  the  square  root  of  the  ratio  between 
the  surface  wind  stress  and  the  air  density.  For  neutral 
stratification,  the  connection  among  the  three  wind  speeds, 
U10,  Ux/2,  and  1/*,  is  the  logarithmic  wind  profile 
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where  Uz  is  the  wind  speed  at  elevation  z,  z0  the  dynamic 
roughness  of  the  ocean  surface,  and  k,  =  0.4  the  von  Karman 
constant.  Equation  (1)  can  also  be  written  as 


e#  =  ct2  gfllAo,  e*  ~  a2  ffttAa,  e**  ~  <r2  g2/ui,  and  so 
on.  The  similarity  relation  of  the  wave  parameters  can 
be  expressed  as 


CT0,5  =  -  In  — ,  (2) 

n  z0 

where  Cz  =  m*/[/2  represents  the  drag  coefficient  referenced 
to  Uz.  Conversion  of  the  growth  functions  from  scaling  with 
U\  o  to  scaling  with  other  wind  velocities  depends  on  a 
reliable  function  of  the  drag  coefficient  of  the  ocean  surface, 
C10.  Kahma  and  Calkoen  [1994]  tested  two  different 
representations  of  C|0  to  convert  the  growth  functions  from 
U[0  to  «*  scaling.  They  reached  the  conclusion  that  the  case 
for  h*  scaling  has  not  been  proved  and  that  if  u*  is  not 
measured,  the  algorithm  chosen  to  estimate  it  may  itself  be  a 
source  of  considerable  error. 

[4]  Based  on  the  consideration  that  surface  waves  are  the 
ocean  surface  roughness  element,  and  that  die  dynamic 
influence  of  surface  waves  decays  exponentially  with  dis¬ 
tance  from  the  air-water  interface,  Hwang  [2004]  processed 
the  drag  coefficient  Cy 2  referenced  to  I/y 2.  The  result 
yields  significant  improvement  in  agreement  among  differ¬ 
ent  data  sets  and  a  similarity  relation  of  Cy 2  dependence  on 
the  dimensionless  frequency  scaled  with  «*  was  developed. 
Subsequently,  the  similarity  property  of  CV2  was  extended 
to  dependence  on  the  dimensionless  frequency  scaled  with 
Uio  [Hwang,  2005a,  2005c],  In  this  paper,  the  wavelength 
scaling  of  drag  coefficient  is  applied  to  die  conversion  of 
wave  growth  functions  from  t/10  to  C/y2  and  u*  scaling.  The 
analytical  expressions  relating  the  three  sets  of  growth 
functions  are  described  in  section  2.  The  analysis  indicates 
that  in  order  to  use  the  similarity  relation  of  Cy2,  the  ratio 
between  Cy2  and  Ul0  (denoted  as  RL/)  needs  to  be  consid¬ 
ered.  A  discussion  of  Rv  is  presented.  Section  3  provides  a 
quantitative  comparison  of  several  different  expressions  of 
the  drag  coefficient.  Section  4  presents  the  wave  growth 
functions  based  on  t/10,  t/y2,  and  u*  scaling.  Data  from 
several  field  measurements  of  fetch-  and  duration-limited 
wave  growtii  studies  and  air-sea  interaction  experiments  are 
compiled  to  compare  with  the  present  conversion  results. 
Comparisons  with  other  proposed  ti*  scaling  functions  as 
well  as  numerical  models  are  also  discussed  in  this  section. 
A  summary  is  given  in  section  5. 


2.  Dimensionless  Wave  Growth  Functions 

[5]  Considering  the  characteristic  parameters  of  the  grav¬ 
ity  wave  field  (wave  height  H„  or  the  variance  of  surface 
displacement  a2,  and  peak  wave  frequency  ljp);  wind  field 
(a  reference  wind  velocity  U),  fetch  x,  duration  t,  and  the 
environment  (water  depth  d,  and  gravity  g),  the  following 
dimensionless  parameters  can  be  established 


a2g2 


Uu0 


"  6  ,  ,/  _  v<  _  Sx  A  _  5‘  J,  _  6“  m 

TTT'^  •*  “  772-'  ~T,'d  (3) 


gt 


gd 


U* 


u 2 


u' 


u2' 


where  a  prime  on  a  variable  denotes  its  dimensionless 
representation.  In  the  following  discussions  when  it  is 
necessary  to  distinguish  among  normalizations  with  t/i0, 
U\j2,  and  «*,  the  dimensionless  variables  are  differ¬ 
entiated  by  subscripts  #,  *,  and  **,  respectively;  that  is, 


J  =f2{xJ,t’,d').  (4) 

[r>]  If  the  local  water  depth  is  sufficiently  deep  and  the 
duration  of  wind  sufficiently  long,  the  fetch-limited  growth 
of  wind -generated  waves  can  be  expressed  simply  by 

e'=f{x),  w'=f2{x!).  (5) 

Extensive  efforts  have  been  devoted  to  the  establishment  of 
the  fetch-limited  growth  functions  f  and  f2  (e.g.,  reviews  by 
Kahma  and  Calkoen  [1992,  1994]  and  Young  [1999]). 
These  efforts  led  to  the  conclusion  that  f\  and  f2  can  be 
represented  by  a  power-law  function  for  a  wide  range  of  the 
dimensionless  fetch,  approximately  102  <  x#<  10  , 

e#=dx"#,  U >#=54.  (6) 

A  summery  of  the  coefficients  A,  a,  B,  and  h  from  several 
different  proposed  functions  is  given  in  Appendix  A.  For.r#  > 
104,  the  rate  of  wave  development  obviously  slows  down 
and  (6)  overestimates  the  dimensionless  wave  variance  or 
wave  period.  Donelan  et  al.  [1992]  developed  a  differential 
expression  of  wave  growth  and  produced  analytical  solutions 
of  the  growth  functions  in  terms  of  x#(w#)  and  x#{e#), 

#  Vl  -  5.5414 ef  2) 

-  2.2690  x  10s  (l  +  (7a) 

x#  =  4.0946  x  104ln[ - — | 

—  0.8302 J 

-  3.3992  x  104(w#  +  0.415l)w!.  (7b) 

As  commented  by  Young  [  1999],  the  transcendental  nature  of 
equations  (7 a)  and  (7b)  makes  solving  e #  and  w#  for  given  x#  a 
cumbersome  process.  Hwang  and  Wang  [2004]  developed  a 
higher-order  data-fitting  technique  to  represent  wave  growth 
in  the  conventional  power-law  functions  (equation  (6))  but 
with  the  wave  development  rate  varying  with  the  dimension¬ 
less  fetch,  duration,  or  wave  frequency.  The  procedure  was 
applied  to  five  different  field  data  sets  obtained  in  steady  wind 
conditions  [ Burling ,  1959;  Hasselmann  etal.,  1973;  Donelan 
etal.,  1985 \ Dobson  etal.,  1989;  Babanin  andSoloviev,  1998] 
(hereafter  referred  to  as  the  BHDDB  data  set;  further 
description  of  the  data  set  is  given  by  Hwang  and  Wang 
[2004]  and  in  Appendix  A).  The  coefficients  for  the  first-  and 
second-order  fitted  growth  functions  are  listed  in  Appendix  A 
(equations  (Al  2)  and  (Al  3)).  With  the  power-law  representa¬ 
tion,  it  becomes  much  easier  to  transform  the  growth 
functions  into  different  dependent  variables.  For  example, 
Hwang  and  Wang  [2004]  converted  the  fetch-limited  growth 
functions  into  duration-limited  growth  functions 

e#=Ptp#,  w#  =  Qtq#.  (8) 
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Figure  1.  The  velocity  ratio  R(r  =  U\^HJ\0,  plotted  as  a  function  of  (a)  to**,  and  (b)  u>#  for  several 
different  wind  speeds.  The  majority  of  reported  10**  and  u>#  falls  between  the  two  vertical  dashed  curves. 


The  coefficients  P,  p,  Q,  and  q  are  simple  algebraic  functions 
of  A,  a,  B,  and  b. 


P  =  A 


[/?,(/>+ 1) 


6+1  « 
,  P-  — 


r  l 
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BbRc(b+  1) 


i+r 


9  = 


b+\' 

1 

b 


(9) 


6+1 


[8]  It  is  difficult  to  obtain  a  consistent  parameterization  of 
the  ocean  surface  drag  coefficient  in  terms  of  Cjo-  Hwang 
[2004,  2005a,  2005b,  2005c]  showed  that  the  similarity 
relation  of  ocean  surface  drag  coefficient  exists  in  wave¬ 
length  scaling,  Cyj2.  To  make  use  of  the  Cy2  similarity  to 
convert  Ui0  scaling  to  m*  scaling  of  the  wave  growth 
functions,  the  following  equalities  can  be  used 


where  Rc  «  0.4  [Yefimov  and  Babanin,  1991]  is  the  ratio 
between  the  effective  group  and  phase  velocities  of  the  wave 
component  at  the  spectral  peak.  The  converted  duration- 
limited  growth  functions  compare  very  well  with  the  limited 
amount  of  data  on  duration-limited  wave  growth  collected 
from  the  ocean  (see  Figure  5  of  Hwang  and  Wang  [2004]). 
Because  reliable  duration-limited  growth  data  are  difficult  to 
acquire  from  the  field,  making  use  of  fetch-limited  growth  data 
to  study  the  temporal  growth  of  surface  waves  is  very  useful. 

[7]  When  the  reference  wind  speed  is  changed  from  U\o 
1°  U\a  or  «*,  the  dimensionless  parameters  are  related  to 
each  other  by  the  drag  coefficient,  Cj0  or  Cy2  [e.g.,  Komen 
et  ai,  1994] 

e**  =  C]02e#,  io#*  =  x**  =  C 10'x#,  f**  =  C](j 

~  =  CyjUJ*  t  ■*’**  =  ■  A*  =  ^\/2^ 

(10) 

Substituting  equation  (10)  to  (6),  the  growth  functions  in 
terms  of  «*  are 


e*  —  Ri*eih  1  x*  —R(fx# >  {*  (12) 


where 


Ru 


U\/2  In  n  —  In  [kpZp) 

U\o  ~  (  “W10  \  ’ 


(13) 


and  kp  is  the  wavenumber  of  the  spectral  peak  component. 
The  dimensionless  roughness  is  related  to  Cy2  by  the 
logarithmic  wind  profile,  equation  (2), 


kpZo  =  it  exp 


(14) 


Ru  can  be  solved  iteratively  for  given  w**  and  U]0.  From 
numerical  experiment,  with  the  initial  guess  of  Ruo  =  1,  a 
relative  error  of  1  percent  is  achieved  within  five  iterations 
[Hwang,  2005c],  Combining  equations  (6),  (10),  and  (13), 
then 


£**  —  cal0  Ax% 


J**  —  C10 


■6+0.5  n..b 

D**  >Jc  • 


(H) 


_  r-a— 2  n2a- 

UX/2  KU 


4Ax%, 


.  _  /^6+0.5  r>26+l  n  b 

u**  —  ^\/2  a£/  d a** 


(15) 


The  duration-limited  growth  functions  in  terms  of  u*  can  be 
derived  from  equation  (11)  in  a  similar  procedure  as  that  of 
deriving  equation  (8)  from  equation  (6). 


Again,  derivation  of  the  duration-limited  growth  function 
from  equation  (15)  is  similar  to  that  of  deriving  equation  (8) 
from  equation  (6). 
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Figure  2.  C10  calculated  using  the  similarity  relation  of  Cxj2  from  wavelength  scaling  [Hwang,  2005c], 
The  DMAJ  data  are  superimposed  for  comparison, 


[9]  Equation  (15)  indicates  that  two  factors,  Cyj2  and  Rv, 
are  needed  to  convert  the  growth  functions  from  Ui0  to  «* 
scaling.  Presently,  the  effort  of  conversion  focuses  on 
equation  (11),  however,  the  uncertainty  in  C10  parameteri¬ 
zation  has  limited  the  progress.  The  issue  of  the  parameter¬ 
ization  of  ocean  surface  drag  coefficient  has  been 
extensively  discussed  in  die  literature  and  a  brief  review 
is  given  in  section  3.  In  section  4,  the  conversion  using 
equation  (15)  will  be  discussed.  Here  the  factor  Ry  is 
examined.  Figure  1  plots  Ry  as  a  function  of  w.*  and  w# 
with  U\  0  as  a  parameter.  Deviation  of  Rv  from  unity  is  a 
source  of  error  in  converting  wave  growth  functions  from 
C/10  to  u*  or  Uyj2  scaling.  As  shown  in  Figure  1,  Rv  varies 
with  dimensionless  frequency  and  wind  speed  in  a  system¬ 
atic  but  complex  fashion.  Interestingly,  as  waves  become 
more  well-developed,  Ry  approaches  asymptotically  to  a 
value  about  1.10  for  U\0  =  5  m/s  and  about  1.35  for  C/I0  = 
20  m/s,  with  a  mean  value  close  to  1 .25  over  the  whole 
range  of  wind  speeds  used  in  the  computation.  Empirically, 
it  has  been  suggested  that  at  full  development,  the  phase 
speed  of  the  wave  component  at  the  spectral  peak,  cp, 
travels  faster  than  the  reference  wind  speed,  C/|0,  with 
Cp/U] 0  ~  1.25  [Pierson  and Moskowitz,  1964],  Using  Ux/2  as 
the  reference  wind  speed,  the  ratio  cpIU\i2  would  have 
been  close  to  unity.  This  provides  an  independent  support 
for  lJ',,/2  serving  as  the  free-stream  velocity,  U^. 


3.  Drag  Coefficient  of  Wind-Generated  Seas 

[10]  Kahma  and  Calkoen  [1992,  1994]  tested  two  differ¬ 
ent  representations  of  C10  to  convert  the  growth  functions 
from  t/,o  scaling  to  «*  scaling.  The  first  one  is  a  modified 


tVu  [1980]  formula  expressing  C]0  as  a  linear  function  of 

C,o, 

C,o  =  8  x  10-4  +  6.5  x  10~5  x  max{{/io,7.5m/s}.  (16) 

The  second  one  is  a  wave-dependent  z0  proposed  by 
Donelan  [1990], 

2  =  5.53x,0-<(^)“  07) 

Figure  2  compares  the  drag  coefficient  computed  from 
equations  (16)  and  (17)  with  the  results  obtained  from  direct 
wind  stress  measurements  in  the  field  under  fetch-limited 
wave  conditions.  The  field  data  set  represents  the  combined 
results  of  four  different  experiments  [ Donelan ,  1979;  Merzi 
and  Graf,  1985;  Audi l  and  Donelan,  1996;  Janssen,  1997] 
(hereafter  referred  to  as  the  DMAJ  data  set;  the  experi¬ 
mental  conditions  were  summarized  by  Hwang  [2004]),  that 
covers  a  wide  range  of  the  wave  development  condition, 
0.0235  <  Wpi<*/g  <  0.237  and  0.0263  <  udcp  <  0.237.  As 
shown  in  Figure  2,  the  Wu  [1980]  expression  of  C,0  is  only 
suitable  for  describing  the  ocean  surface  drag  condition  at 
more  mature  sea  state  (for  w#  in  the  neighborhood  of  unity). 
C, 0  computed  from  Donelan  [1990]  roughness  expression, 
equation  (17),  shows  dependence  on  both  C/,0  and  w#  but 
the  range  of  variation  on  in#  is  smaller  than  that  observed  in 
the  field  data  by  about  a  factor-of-two  (Figure  2).  Hwang 
[2004,  2005a,  2005b,  2005c]  suggested  that  the  difficulty  in 
finding  the  similarity  properties  of  the  ocean-surface  drag 
coefficient  can  be  attributed  to  the  choice  of  the  arbitrary 
10-m  as  the  length  scale  for  wind-speed  reference.  When 
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Figure  3.  Fetch-  and  duration-limited  wave  growtli  data  represented  as  with  (a)  U]0,  (b)  £/v 2>  and 
(c)  u*  scaling. 


processed  with  wavelength  scaling,  a  strong  similarity 
behavior  exists  in  wind-stress  measurements  represented 
either  by  the  drag  coefficient  or  the  dynamic  roughness, 

~~J  i  with  Ac  =  1.220  x  10~2  and  ac  =  0.704. 

(18) 

Using  equation  (18),  the  dimensionless  roughness  expressed 
as  /r;jz0  can  be  easily  computed  with  the  application  of  the 
logarithmic  wind  speed  profile  (equation  (14)).  The 
dimensionless  roughness  can  also  be  expressed  as  z0At 
with  the  additional  application  of  the  fetch-limited  growth 
function,  e#(w#).  Furthermore,  C10  (U10,w#)  =  C^2  Ru  can 
be  derived  from  equations  (18)  and  (13)  [Hwang,  2005c]. 
The  result  is  in  good  agreement  with  field  observations 
(Figure  2).  For  practical  applications,  the  following 
parameterization  function  is  developed  from  the  DMAJ 
data  set  [Hwang,  2005a,  2005c] 

C\/2  =/lu,(U1/'g10')  'With  ,4|o  =  1.289  x  10-3  and  H|0  =0.815. 

(19) 


4.  Comparison  with  Field  Data  and  Model 
Results 

[n]  Most  experiments  on  fetch-limited  wave  growth 
studies,  such  as  BHDDB,  do  not  report  direct  wind  stress 
measurements  so  they  cannot  provide  an  unambiguous 
assessment  on  the  growth  functions  using  different  scaling 
wind  velocities.  Direct  wind  stress  measurements  are  usu¬ 


ally  performed  in  air-sea  interaction  experiments,  especially 
for  the  parameterization  of  ocean  surface  drag  coefficient. 
The  DMAJ  data  set  was  obtained  under  fetch-limited 
growth  conditions  with  sufficient  details  reported  to  facili¬ 
tate  the  computation  of  Uyj2-  Figure  3  shows  the  comparison 
of  these  two  large  data  sets.  Because  many  of  the  experi¬ 
ments  on  air-sea  interaction  do  not  report  wind  fetch,  the 
comparison  is  shown  in  e'(w').  In  addition  to  the  two 
combined  fetch-limited  data  sets,  measurements  from  a 
recent  duration-limited  wave  growth  study  [Hwang  and 
Wang,  2004]  are  also  superimposed.  Considering  the  dis¬ 
parate  conditions  under  which  these  data  were  collected,  the 
varieties  of  instruments  used,  and  tire  differences  in  pro¬ 
cessing  procedures  by  different  research  groups,  the  general 
agreement  of  the  results  shown  in  Figure  3a  is  impressive. 
The  curve  plotted  in  the  panel  used  for  visual  reference  is 
equation  (A  13),  which  is  derived  from  the  second-order 
fitting  procedure  applied  to  the  BHDDB  data  set  [Hwang 
and  Wang,  2004].  In  Figures  3b  and  3c,  equation  (15)  is 
used  to  convert  the  wave  growth  functions  from  U10  to  Uyj2 
and  u*  scaling.  As  mentioned  earlier,  to  use  equation  (15), 
C\j2  and  Ru  need  to  be  calculated.  For  the  DMAJ  data  set, 
direct  measurements  of  u*  and  kp  were  reported  and  the 
computation  of  Cyj2  and  Ru  is  straightforward  [Hwang, 
2004,  2005c],  For  the  duration-limited  data  set  of  Hwang 
and  Wang  [2004],  kp  is  available,  but  the  fast  evolution  of 
the  wave  field  requires  high  temporal  resolution,  on  the 
order  of  164-s  duration  for  each  wave  spectral  computation. 
Such  duration  is  too  short  for  reliable  wind  stress  processing 
and  bulk  formulae  are  used  to  obtain  u*.  Here  computed 
by  equation  (19)  is  chosen.  More  details  on  «*  computation 
using  bulk  formulae  are  given  in  Appendix  B.  For  the 
BHDDB  data  set  (digitized  from  dimensionless  plots  of 
e#(xfj)  or  e#( w#)  and  w#(x#)  in  the  original  papers),  simulta- 
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Figure  4.  (a)  Fetch-,  and  (b)  duration-limited  growth  functions  in  u*  scaling,  and  comparison  with 
available  field  data. 


neous  w#  and  Ul0  are  not  available  to  compute  RL,  (Figure  1).  the  DMAJ  data  set  with  direct  u*  measurements,  and  the 
The  curves  corresponding  to  Uu)  =  7  and  20  m/s  are  applied  Hwang  and  Wang  [2004]  data  set  with  u*  calculated  by 
to  equation  (A13)  to  represent  the  possible  range  of  the  equation  (19). 

average  of  the  BHDDB  data  set.  The  growth  functions  [12]  The  growth  functions  e**(x**),  w**(x**),  e**(t**), 
scaling  with  UW2  311(1  »*  and  converted  from  those  with  Uw  and  w**(f**)  are  shown  in  Figure  4.  A  subset  of  the  DMAJ 
scaling  are  in  very  good  agreement  with  the  results  based  on  data  set,  Donelan  [1979]  and  Merzi  and  Graf[  1985],  listed 


Figure  5.  Comparison  of  the  wave  growth  functions  with  several  other  different  expressions  and 
numerical  model  results,  (a)  C/10  scaling,  fetch-limited  growth;  (b)  w*  scaling,  fetch-limited  growth;  and 
(c)  «*  scaling,  duration-limited  growth. 
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wind  fetch  so  they  can  be  compared  with  the  analytical 
computation  of  the  fetch-limited  wave  growth  functions 
(Figure  4a).  The  duration-limited  measurements  of  Hwang 
and  Wang  [2004]  are  shown  in  Figure  4b  to  compare  with 
the  analytical  computation  of  the  duration-limited  wave 
growth  functions.  The  agreement  appears  to  be  very  good 
considering  the  large  data  scatter  in  the  measurements. 

[13]  Figure  5  presents  the  wave  growth  functions  scaled  by 
C/10  and  «*.  In  panel  (a),  several  proposed  growth  functions 
of  e#(x#)  and  u#(x#),  as  summarized  in  Appendix  A,  are 
plotted  together  [. Hasselmann  et  al.,  1973;  Donelan  et  al., 
1992;  Kahma  and  Calkoen,  1994;  Hwang  and  Wang, 
2004].  In  general,  all  different  expressions  of  the  growth 
functions  show  very  good  agreement  in  the  middle  range  of 
the  dimensionless  fetch  about  x#  =  103.  Deviation  of  the 
simple  power-law  growth  functions  from  observations  is 
found  at  earlier  and  later  stages  of  wave  development.  The 
slowing  down  of  wave  growth  as  dimensionless  fetch 
increases  is  reproduced  reasonably  well  by  numerical  wave 
models  [e.g.,  Janssen  et  al.,  1994;  Janssen,  2004].  The 
disadvantage  of  using  Ui0  as  the  scaling  wind  speed  is 
demonstrated  clearly  in  numerical  experiments.  An  example 
is  given  in  Figure  3.22  of  Janssen  et  al.  [1994],  which  is 
reproduced  in  Figure  5a,  showing  that  the  computed  di¬ 
mensionless  function  ui#(x#)  displays  additional  dependence 
on  the  wind  speed.  Scaling  with  u*  is  considered  to  be  a 
solution  to  the  problem  associated  with  C/,0  scaling.  Be¬ 
cause  the  uncertainty  in  finding  a  suitable  C,0  parameter¬ 
ization,  it  remains  a  difficult  task  establishing  growth 
functions  using  u*  scaling.  The  results  reported  by  Kahma 
and  Calkoen  [1994],  as  described  in  section  3,  are  shown  in 
Figure  5b.  Also  plotted  in  this  panel  are  the  numerical 
model  results  reported  by  Janssen  et  al.  [1994,  Figures  3.25 
and  3.26]  and  Janssen  [2004,  Figure  5.10]  and  the  conver¬ 
sion  using  the  Cy2  similarity  relation  described  in  section  2. 
Figure  5c  shows  the  result  on  duration-limited  growth 
functions.  The  numerical  results  of  duration-limited  growth 
given  by  Janssen  [2004,  Figure  5.7]  are  superimposed. 
Discrepancies  between  numerical  simulations  and  the  em¬ 
pirical  functions  of  the  present  analysis  can  be  attributed  to 
several  factors.  First,  differences  in  the  data  sets  used  for 
constructing  the  growth  functions  that  the  wave  model  was 
tuned  to.  Second,  the  drag  coefficient  used  in  converting 
from  C/,o  to  «*  scaling;  more  detail  discussion  of  the  issue 
is  given  in  Appendix  B.  Briefly,  e**  is  proportional  to  «*4, 
which  in  turn  is  proportional  to  C72.  Therefore,  a  factor-of- 
two  underestimation  of  C7  (not  uncommon  for  young  wave 
fields,  see  Appendix  B)  would  overestimate  e**  by  a  factor- 
of-four.  This  is  not  only  a  problem  for  duration-limited 
growth  study,  it  also  occurs  in  fetch-limited  growth,  as 
illustrated  by  the  large  difference  between  the  two  numer¬ 
ical  model  curves  computed  for  8  and  18  m/s  winds 
[Janssen  et  al.,  1994]  shown  in  Figure  5a.  Third,  the 
incorporation  of  Ry  in  the  conversion  method  also  plays  a 
minor  role.  A  fourth  factor  contributing  to  the  discrepancy 
between  observed  and  modeled  duration-limited  growth 
curves  is  the  initial  condition  (of  peak  spectral  frequency 
at  0.34  Hz)  and  the  resolved  frequency  range  (maximum 
frequency  0.7  Hz)  of  model  runs  used  in  Janssen  [2004], 
Janssen  (personal  communication,  2005)  performed  the 
same  numerical  experiment  of  duration-limited  growth  with 
a  higher  initial  peak  frequency  (0.9  Hz)  and  the  frequency 


range  of  computation  expanded  to  2  Hz.  The  trend  of  the 
wave  growth  (shown  as  circles  in  Figure  5c)  is  in  much 
better  agreement  with  the  empirical  curve  derived  from  the 
present  analysis. 

[14]  Because  fetch-  and  duration-limited  growth  data  with 
n*  measurements  are  rare,  conversion  of  e#(x#)  and  oj#(x#)  to 
e**(x**),  w**(x**),  e**(t**),  and  u>**(/**)  remains  a  nec¬ 
essary  task  to  enlarge  the  field  database  for  model  validation 
and  verification.  The  method  described  in  this  paper  makes 
use  of  the  similarity  relation  of  the  drag  coefficient  based  on 
wavelength  scaling  and  considers  the  factor  t/y2/t/|o- 

5.  Summary 

[15]  While  C/,o  is  used  extensively  as  the  reference  wind 
speed  in  the  research  of  wind-wave  growth  and  air-sea 
interaction,  the  choice  of  10-m  as  the  reference  elevation 
of  wind  speed  measurement  is  mainly  due  to  practical 
consideration  rather  than  the  dynamic  significance  of  the 
10-m  elevation  in  the  marine  boundary  layer.  Searching  for 
an  alternative  scaling  velocity  has  been  a  continuous  effort 
over  the  last  few  decades.  It  is  believed  that  either  [4/2  or  u* 
is  more  preferable  than  Ul0  to  serve  as  the  scaling  wind 
speed.  The  former  is  a  reasonable  substitute  for  the  free- 
stream  velocity  used  in  the  dimensionless  analysis  of  wave 
growth  ftmctions  (section  2),  and  the  latter  represents  the 
wind  stress  applied  at  the  ocean  surface.  Converting  from 
f/,o  to  U\a  or  11  *  scaling  requires  an  accurate  prescription  of 
the  drag  coefficient  of  the  ocean  surface.  Although  there 
have  been  many  fonnulae  proposed  for  C,0,  it  is  difficult  to 
express  its  complex  dependence  on  Ul0  and  in  an 
analytical  form  amenable  to  computation.  Recently,  Hwang 
[2004,  2005a,  2005b,  2005c]  showed  that  the  parameteri¬ 
zation  of  ocean  surface  drag  coefficient  simplifies  consid¬ 
erably  when  [/V2  is  used  as  the  velocity  scale,  and  similarity 
relation  exists  in  both  the  drag  coefficient  and  the  dynamic 
roughness.  The  similarity  relation  of  the  drag  coefficient 
based  on  wavelength  scaling  is  used  in  this  paper  to  convert 
the  wave  growth  functions  from  C/,0  to  64/2  and  w*  scaling. 
The  results  are  in  good  agreement  with  the  limited  field  data 
that  provide  w*  measurements  and  fetch  or  duration  infor¬ 
mation  to  form  dimensionless  parameters  scaled  with  w*. 

Appendix  A:  A  Summary  of  the  Fetch-Limited 
Growth  Functions 

[te]  The  following  fetch-limited  growth  functions  are 
mostly  compiled  by  Kahma  and  Calkoen  [1994]  and  Young 
[1999], 

[17]  Sverdrup,  Munk  and  Bretschneider  empirical  func¬ 
tions  (SMB)  [CERC,  1977] 

e#  =  5.0  x  10“3  tanh2  (o.0125xi42Y  u)#  = - - v. 

V  *  >  *  tanh(0.077425) 

(Al) 

(Note:  combined  laboratory  and  field  data  sources.) 

[is]  Pierson-Moskowitz  limit  of  fully-developed  seas 
[Pierson  and  Moskowitz,  1964] 

e#  =  3.64  x  10~3,  u>#  =  0.82.  (A2) 
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Table  Al.  Coefficients  of  the  Fetch-Growth  Similarity  Laws 


l’d  =  A  Xtf  and  u#  —  B  x#h 

Source 

A 

a 

B 

b 

1. 

SMB  [ CERC ,  1977]“ 

7.82  x  10"7 

0.84 

10.82 

-0.25 

2. 

JONSWAP  [Hasselmcmn  et  al.,  1973]“ 

1.60  x  10“7 

1.00 

21.98 

-0.33 

3. 

Bothnian  Sea  (unstable)  [Kahma,  1981] 

3.60  x  10"7 

1.00 

19.97 

-0.33 

4. 

Lake  Ontario  [Donelan  et  al. ,  1985] 

8.42  x  10"7 

0.76 

11.62 

-0.23 

5. 

North  Atlantic  [Dobson  et  al.,  1989] 

1.27  x  10“6 

0.75 

10.68 

-0.24 

6. 

Lake  St.  Clair  [Donelan  et  al. ,  1992]b 

2.60  x  10“7 

0.95 

17.59 

-0.30 

7. 

Composite  stable  [Kahma  and  Calkoen,  1994] 

9.30  x  10~7 

0.77 

12.00 

-0.24 

a. 

Composite  unstable  f Kahma  and  Calkoen,  1994] 

5.40  x  10-7 

0.94 

14.00 

-0.28 

9. 

Composite  mixed  [Kahma  and  Calkoen ,  1994] 

5.20  x  10~7 

0.90 

13.70 

-0.27 

10. 

Average  [Young,  1999] 

7.50  x  10~7 

0.80 

12.56 

-0.25 

11. 

Fetch-dependent  growth  rate  (equations  (A  12)  and  (A  13))  \ Hwang  and  Wang,  2004] 

6.19  x  10"7 

0.81 

11.86 

-0.24 

“SMB  and  JONSWAP  combine  field  and  laboratory  data. 

bThe  asymptotic  form  of  the  original  implicit  functions  of  dimensionless  energy  and  frequency  for  the  range  1 00  <  x*  <  3000. 


[19]  JONSWAP  fetch-limited  wave  evolution  [Hasselmann 
et  al,  1973] 

e#  =  1.6  x  l(T7x#,  uj#  =  21.98  x/33.  (A3) 

(Note:  combined  laboratory  and  field  data  sources.) 

[20]  Bothnian  Sea  [Kahma,  1981] 

e#  =  3.60  x  1(T7x#,  w#  =  19.97a/'33.  (A4) 

(Note:  strongly  unstable  conditions.) 

[21]  Lake  Ontario  [Donelan  et  al.,  1985] 

eit.  =  2.74  x  10-3w/'3,  w#  =  11.6  a/23.  (A5a,A5b) 

Combining  (A5a)  and  (A5b), 

e#  =  8.415  x  KT’x/6.  (A5c) 

Nortli  Atlantic  open  ocean  [Dobson  et  al,  1989] 

e#  =  1 .27  x  l(T6x/5,w#  =  10.68  a/24.  (A6) 


(Note:  the  analysis  took  into  account  the  coastal  effect  on 
the  wind  field  and  used  the  integrated  wind  speed  up-fetch 
as  the  scaling  wind  velocity.) 

[22]  Lake  St.  Clair  [. Donelan  et  al,  1992] 


a#  =4.0946  x  104  In 


1  -  5.5414  e'f  2 


1/32^1/3.2 
'#  ’ 


2.269  x  105(l  +  2.7707  e/'2)e, 


(  UJ*  ^ 

—  0.829/ 


a#  =4.0946  x  104hi 


(ur#  —  0.829 ) 
8.616  x  10-4(w#  +  0.414) 


4 


(A7a) 


(A7b) 


[23]  Composite  field  data  [Kahma  and  Calkoen,  1992, 
1994], 

[24]  The  Bothnian  Sea  data  set  highlights  the  significant 
influence  of  the  stability  effect.  Kahma  and  Calkoen  [1992, 
1994]  suggested  the  following  stability  adjustment  based  on 
reanalysis  of  several  field  data  sets. 

[25]  Stable  stratification: 


e#  =  9.3  x  10“7.v|77,  u>#  =  12  a/24.  (A8) 

Unstable  stratification: 


e#  =  5.4  x  ]0~7x/4,  w#  =  14  a/28.  (A9) 

Combined  field  data: 

e#  =  5.2  x  10-7a°/,  w#  =  13.7  a/'27.  (A10) 


Average  and  bounds  [Young,  1999] 


e#  =  max 


{ 


(7.5  ±  2.0)  x  IO-7*#8 
(3.6  ±  0.9)  x  10-3  ’ 


(12.56  ±  1.88)x/25 
(0.82  ±0.13) 


(All) 


Fetch-dependent  growth  rate  [Hwang  and  Wang,  2004], 
[26]  Hwang  and  Wang  [2004]  included  the  data  sets 
reported  by  Burling  [1959]  and  Babanin  and  Soloviev 
[1998]  with  other  measurements  (BHDDB)  to  derive  the 
fetch  laws 


e#  =  6.1910  x  l(T7x/106,  w#  =  11.86  a/2368.  (A12) 

They  also  developed  a  higher-order  fitting  technique  to 
describe  wave  growth  in  power-law  functions  with  variable 
proportionality  coefficient  and  exponent.  To  the  second 
order. 


(Note:  eif  and  are  implicit  functions  of  xtl.  The  additional 
digits  beyond  second  decimal  place  in  (A7a)  and  (A7b)  not 
available  in  the  original  paper,  were  supplied  by  M. 
Donelan  (personal  communication).) 


e#  =  42xg,  U)#=B2*#, 

A2  =  exp(  «o)a/2  lnj'# ,  a2  =  00  +  2a2lnx# , 
B2  =  exp(30)x/2  ,  b2  =  (3,  +  2p2lnx#, 


(A13) 
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Figure  Bl.  (a)  C10  of  Hwang  and  Wang  [2004]  data  set  calculated  using  the  similarity  relation  of  Cy/2 
(equation  (19)),  wind-dependent  Cw  [Wu,  1980],  and  the  LKB  bulk  parameterization  [Liu  et  al.,  1979], 
The  DMAJ  data  are  superimposed  for  comparison,  (b)  The  impact  on  the  growth  curves  by  using 
different  u*. 


with  a0  =  3.0377,  a,  =  -0.3990,  a2  =  0.0110,  (30  = 
-17.6158,  (3,  =  1.7645,  and  P2  =  -0.0647.  Table  Al 
presents  a  summary  of  A,  a,  B,  and  b  of  various  growth 
functions  described  above. 

Appendix  B:  Computation  of  h*  of  Hwang  and 
Wang  [2004]  Data  Set 

[27]  The  duration-limited  wave  growth  data  reported  by 
Hwang  and  Wang  [2004]  were  collected  in  the  first  two 
hours  of  wave  development  after  the  beginning  of  a  steady 
wind  event  in  a  sheltered  bay.  Because  the  wave  develop¬ 
ment  is  very  fast  at  such  an  early  stage,  the  duration  of  data 
segments  for  wave  spectral  processing  is  about  1 64  s.  This 
duration  is  too  short  for  reliable  wind  stress  computation 
and  parameterization  equations  need  to  be  used  to  calculate 
«*.  Figure  Bla  shows  the  results  obtained  from  the  LKB 
bulk  parameterization  [Liu  et  al.,  1979],  wind-speed  depen¬ 
dent  C)0  by  Wu  [1980],  and  Cv2  similarity  function, 
equation  (19).  The  trend  of  C10(t/io>  w#)  derived  from 
equation  (19)  is  consistent  with  the  results  displayed  by 
the  direct  wind  stress  measurements  of  the  DMAJ  data  set. 
In  contrast,  the  LKB  and  Wu  parameterizations  seem  to  be 
unsuitable  for  application  to  very  young  seas.  Because  e** 
is  proportional  to  w* 4  and  f**  is  proportional  to  h*  ,  when  a 
lower  value  is  used  in  the  scaling,  the  growth  curve  shifts 
upward  and  right-ward,  but  the  upward  shift  is  four  times  as 
much  as  the  right-ward  shift  in  log-log  scales  (Figure  Bib). 
In  Figure  Bib,  the  numerical  model  results  of  e**(f**) 
reported  by  Janssen  [2004,  Figure  5.7]  and  the  rerun 
(Janssen,  personal  communication,  2005)  as  described  in 
the  text  are  also  shown. 
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